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Combined Theory of Reflectance 
and Emittance Spectroscopy 


At visible to infrared wavelengths, interaction of electromagnetic radiation with a 
particulate medium involves a combination of reflection, absorption, scattering, and 
emission processes, each dependent on the geometry of the interaction, the physical 
character of the surface, and the wavelength of radiation being analyzed. Approximate 
analytic expressions provide a useful framework to evaluate individual components and 
their potential effects on spectral features. 

Chapter contributed by Bruce Hapke. 


/ In-this-chapter the theory in which either or both 
reflected sunlight and thermally emitted radiation con- 
tribute to the power received by a detector viewing 
a particulate medium, such as a powder in the 
laboratory or a planetary regolith, is considered theo- 
retically. This theory is of considerable interest for the 
interpretation of data from field or spacecraft instru- 
ments that are sensitive to the near-infrared region 
of the spectrum, such as NIMS (near-infrared mapping 
spectrometer) and VIMS (visual and infrared mapping 
spectrometer), as well as thermal infrared detectors. 
Exper^mental .andiObs^tyatjonal techniques and re- 
sults are discussed extensively in several other chap- 
ters -in -this-book. 

There exists a large body of literature that treats 
radiation transfer in planetary and stellar atmospheres 
(eg.. Chandrasekhar, i960; Goody, 1964;Sobe/oy 1975; 
Van de Hulst, 1980), but it has not been widely recog- 
nized until recently that models based on the equation 
of radiative transfer also apply to powdered materials 
in which the particles are closely spaced. However, 
the exact theories require the solution of nonlinear 
integral equations. Since in most cases the absolute 


radiance received and the particle scattering proper- 
ties of the surface being investigated are poorly known, 
exact theories are not necessarily more useful, and are 
much less convenient, than approximate analytical 
solutions. 

Theoretical treatments of reflectance or emittance 
from a particulate medium have been published by 
Vincent and Hunt (1968), Conel (1969), Emslie and 
Aronson (1973), Lumme and Bowel! (1981), and Hapke 
(1981, 1986). However, the situation when both re- 
flectance and emittance contribute to the received 
signal has not been discussed. Also, most theories of 
emittance by powders derive only the hemispherically 
integrated radiance, whereas it is the directional 
emittance that is measured from spacecraft. 

In this chapter the author's reflectance model 
(Hapke, 1981, 1 986) is extended to include the effects 
of thermal radiation. The general approach and as- 
sumptions on which that reflectance model is based 
also apply here. The particles of the medium are 
assumed to be larger than the wavelength, irregular 
and randomly oriented and positioned, so that coher- 
ent and near-field effects in the radiation scattered 
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LIST OF SYMBOLS 

AA: surface area 

r: distance from detector 

Aa: detector area 

r dd : bidirectional reflectance 

B: Planck emission function 

r dh : directional-hemispherical reflectance 

B 0 : constant term in emission function 

r hd : hemispherical-directional reflectance 

B,: coefficient in emission function 

ri, h : bihemispherical reflectance 

b: opposition effect function 

S: volume scattering coefficient 

b Q : amplitude of opposition effect 

s: distance 

BRDF: bidirectional reflectance distribution function 

T: temperature 

c: speed of light 

u = Ez: dimensionless distance 

D: particle diameter 

V: volume in medium 

E: volume extinction coefficient 

v: volume of particle 

e: angle of observation 

w: single scattering albedo 

F: volume emission coefficient 

z: altitude 

f: filling factor 

G: volume angular scattering coefficient 

g: phase angle 

y = v 1 - w: albedo factor 

H(w,x): H-function 

e: emissivity 

h: Planck's constant 

e d : directional emissivity 

I: radiance 

£ h : hemispherical emissivity 

i: angle of incidence 

4>: emitted integrated flux 

J: incident irradiance 

<t>: hemispherically integrated flux 

j: subscript indicating type of particle 

0. polar angle 

K: volume absorption coefficient 

A: thermal emission scale height 

k: opposition effect width parameter 

A: wavelength 

L = EA: dimensionless emission scale height 

ft: cose 

n: number of particles per unit volume 

H 0 - cosi 

p: average particle angular scattering function 

a. particle geometric cross section 

Q a : absorption efficiency 

Cl. direction 

Q c : extinction efficiency 

dfl: increment of solid angle 

Q s : scattering efficiency 

Aw: solid angle of detector 


between particles can be neglected. An approximate 
solution for the emergent radiance based on the two- 
stream approximation to the equation of radiative 
transfer will be derived. 

Two important omissions from this model are 
the effects of large-scale surface roughness and the 
neglect of polarization. The changes that roughness 
causes to the reflectance are examined in Hapke 
( 1 984), but the analogous theory for emittance is much 
more difficult and is beyond the scope of this chapter. 
A recent treatment is that of Spencer (1990). The 
neglect of polarization greatly simplifies the deriva- 
tions and can be partially justified on the grounds that 
large, irregular, dielectric particles do not polarize light 
strongly by scattering. 

2, 1 . DERIVATION OF THE 
EMITTED RADIANCE 

The geometry of the model is shown schematically 
in Fig. 2. 1 . The plane surface z = 0 separates an empty 
half-space z>0 from a half-space z<0 filled with 


particles. The medium is illuminated by a collimated 
irradiance J; the direction to the source makes an angle 
i with the upward normal to the surface. The medium 
is observed by a detector whose line of sight, as seen 
from the surface, makes an angle e with the normal. 
The phase angle between the source and detector is 
denoted by g. The detector has sensitive area Aa and 
accepts radiation from within a solid angle Aw. Let 

M = cose (2.1a) 

and 

A/ 0 = cosi ( 2 . 1 b) 

Assume that the medium is made up of different 
types of particles denoted by subscript j. The particles 
may differ because of size, composition, or other 
characteristics. Then the extinction, scattering, and ab- 
sorption coefficients of the medium are respectively 

E SjnjOjQ e j 


(2.2a) 
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Fig. 2.1. The geometry of the model. 


S = SjnjOjQsj (2.2b) 

K = SjnjcrjQaj (2.2c) 


F = XjnjOjCj (2.2e) 

The efficiencies and angular scattering functions of 
isolated particles are treated theoretically in several 
places (e.g.. Van de Hulst, 1957; Kerker, 1969; Bohren 
and Huffman, 1983). Except for Q ej , these quantities 
do not change appreciably when the particles are close 
together (Hapke, 1981). Since Q sj + Q aj = Q ej , then S + 
K = E. In general, E, S, K, G, and F are functions of 
wavelength X. 

Let the volume of particles of type j be Vj. Then the 
filling factor, or fraction of the volume occupied by 
solid matter, is 

f = Sjnp/j (2.20 

and the porosity is 1 - f. Strictly speaking, the densities 
nj in equations (2.2a-e) should actually be replaced 
by the effective density -njln(l-f)/f (Hapke, 1986); 
however, this distinction will be found to have no effect 
on the final form of the equations except for the 
opposition effect, which will be discussed separately 
below. 

Define the average properties of the medium as 
follows 

average single scattering albedo, 
w = S/E=[2 j n j a j Q sj ]/[2 j n j a j Q ej ] ( 3a) 


In these equations nj is the number of particles of type 
j per unit volume, CTj is the rotationally averaged geo- 
metric cross-sectional area of the jth type of particle, 
and Q e j, Q sj , and Q aj denote their extinction, scattering, 
and absorption efficiencies, respectively. 

If a particle is large compared with the wavelength 
X and is isolated, then Q ej = 2. However, when the 
particles are so close together that they touch, the 
diffracted light must be regarded as associated with 
the holes between the particles rather than with the 
particles. Because the particles are close together the 
light diffracted by the holes does not have sufficient 
space to spread appreciably before encountering a 
particle and therefore cannnot be distinguished from 
the incident, unextinguished light. Hence, for particles 
close together, Q ej = 1 . 

Denote the angular scattering functions of each type 
of particle by Pj(g), where p ( (g) is normalized so that 
its integral over all solid angles is equal to 4tt, and 
the emissivity of a particle of type j by t v Then the 
angular scattering coefficient of the medium is 


G(g)= 2 j n j a j Q sj Pj(g) (2.2d) 


and the emission coefficient is 


albedo facto r, 
y = V 1 -w 

average particle angular scattering function, 
P(§) G(g)/S [Sj nj cf j Qsj Pj (g)3 / [2j nj C7j Q S j] 

average particle emissivity, 
e F/E— [SjnjCTj€j]/[SjnjCTjQ e j] 


(2.3b) 


(2.3c) 


(2.3d) 


In general, the parameters w, y, p(g), and e are all 
functions of X. Furthermore, it can readily be shown 
that q = Q aj , so that 


F= K 


(2.3e) 


and 


e = K/E=l-W=y 2 (2.3f) 

Although the detector appears to be examining an 
area AA located on the apparent surface a distance 
R from the detector, it is actually receiving radiation 
from all the particles below the surface within the cone 
Aw. Consider an increment of volume dV = r 2 drAw 
located at depth z below the apparent surface and 
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a distance r from the detector. There are three con- 
tributions to the power received by the detector from 
dV: (1) radiation from the source scattered once by 
the particles in dV into the direction toward the 
detector; (2) radiation thermally emitted by the par- 
ticles in dV toward the detector; and (3) radiation that 
has been emitted or scattered at least once, impinging 
on the particles in dV and being scattered toward the 
detector. After the radiation emerges from dV it is 
exponentially attenuated as it traverses the medium 
on its way to the surface. 

Thus, the power received by the detector can be 
written mathematically in the following form 


AP = 


/ [ 

r = R 


j e -E|z|/*> + l B(T) + 

4iT 7 T 


/ 


G(g') 

i(z,n')-~ dn' 
4rr 


>-E|z|//< 


dV 


Aa 


/ [ 


je-EUkMo wp(g) + 4tB(T) + 


r = R 


wj i(z,n , )p(g , )dn' 


e - e|z|/ji £ r 2(jrAftj 


Aa 

4rrr 2 


4 rr 


(2.4) 


where B(T) = is the Planck blackbody thermal emis- 
sion function 


B(T) = 


2zrhc 2 

A 5 


^ghc/XffT_ | j -1 


(2.5) 


T(z) is the temperature, h is Planck's constant, c is 
the speed of light, « is Boltzmann's constant, I(z,fl) 
is the radiance in the medium at position z traveling 
into direction Cl of radiation that has been emitted 
or scattered at least once, and g' is the angle between 
directions and (1 Note that B(T) is the power emitted 
per unit area, so that the power emitted per unit area 
per unit solid angle is B(T)/tt. Equation (2.4) assumes 
that the emittance of each particle is isotropic, which 
will be true for an ensemble of randomly oriented 
particles even though they may not necessarily scatter 
isotropically. 

Now, the reciprocal of E has units of distance and 
is known as the extinction length. Let 


u = Ez (2.6) 

that is, u is the dimensionless altitude above or below 


the apparent surface measured in units of the ex- 
tinction length. Then, since dr= -dz//z. the radiance 
at the detector emerging from the surface can be 
written 


J(i,e,g) 


AP 

Aa Aw 




Jwp(g)e u/ ''° + 4-y 2 B(T) + 


/ . du 

I(u,n')p(g)dft' e u/> — (2.7) 


where-it has been explicitly recognized that z and u 
are negative over the integral. 

Since the temperature is a function of altitude z, 
the thermal emission function B(T) is also a function 
of z and can be written B(z). Expression (2.7) will be 
evaluated under the condition that B(z) can be ap- 
proximated by 

B(z) = B 0 + B,e' |z|/A (2.8a) 

or 

B(u) = B 0 + B,e‘ |u|/L (2.8b) 

where A is a thermal emission scale height, and L = 
EA. Both B 0 and B, are functions of A. This expression 
is sufficiently general that it can describe a subsurface 
temperature distribution containing a thermal gra- 
dient, yet is simple enough that an analytic solution 
can be obtained using the two-stream method. For 
instance, B(0) = B o + B,, B(-°°) = b 0 , and (<9B/<9z) | 
z = 0 = (<9B/<3T) | z = o(dT/ dz)| z = 0 = -B)/A. 

The first two terms in the integrand of equation (7) 
can be evaluated directly. The two-stream approxi- 
mation to the equation of radiative transfer will be 
used to evaluate the third term. Including a collimated 
source and thermal emission, the equation of radiative 
transfer is 


dI(z,H) 1 r 

= -EI(z,fl) H / I(z,nr)G(g')dft' + 

ds 47 t J 

4 7T 

j e -E|z|/^5® +E B(Z) 

477 7T 

where ds is an increment of length in the direction 
fl. This equation describes the changes dl that occur 
as a ray of radiance 1 traverses a distance ds parallel 
to direction fl in a medium that scatters, absorbs, and 
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emits. The first term on the right describes the energy 
removed from the beam, the second and third terms 
describe the energy added by scattering, and the last 
term describes the contribution of thermal emittance. 

Let 0 be the angle between the upward normal and 
ds (Fig. 2.2), so that ds = dzsecfl. Dividing through by 
E and putting du = Edz, the radiative transfer equation 
can be written in the form 

dl(u,fl) w r 

cose =-I(u,fl)+ — / i(u,ft’)p(flf)dn’ + 

du 47 r J 

An 

w y 2 

) — p(g)e u/ ^ + — B(u) (2.9)' 

477 7r 


seen by the following. Suppose p(g) scatters light much 
more strongly in the backward direction than in the 
forward. Then I will be stronger in the upward direction 
than if p(g) were isotropic; however, this excess will 
be preferentially scattered downward when I interacts 
with p(g), so that the effects of the anisotropy on the 
third term will approximately cancel each other. 

Hence, in the first approximation, provided the de- 
partures from isotropy in p(g) are not too great, in 
equation (2.7) we may put 

I(unr)p(gf)dft « J Ij(ufl')dn' (2.10) 

4;r 4 t r 


Now, the third term in the integrand of equation (2.7) 
refers to the multiply scattered component of the 
radiance field in the medium. This term will be much 
less sensitive to p(g) than the first term, as can be 


Z 



Fig. 2.2. The equation of radiative transfer describes the 
changes dl that occur in the beam of radiant energy as it 
traverses the cylinder of length ds. 


where 1| (u,H) is the solution of equation (2.9) for iso- 
tropic scattering. If the scatterers are isotropic, p(g) = 
1, and equation (2.9) is 


dMuft) w r 

cose- — - — =-ij(utt) + — / Munidn' + 
du 4ir J 

An 


w . y 

J — e u/Mo + — B(u) (2.11) 

477 77 


Using the two-stream method for solving equa- 
tion (2.1 1), the equation is integrated separately over 
the upward-going and downward-going hemispheres, 
and I; is replaced by its appropriate average in the 
integrals. Thus, let 4>u and <f> 0 be the integrated power 
per unit area traveling into the upward and downward 
directions, respectively 


= J Mun'JdfT = J 


</>u(u) = j Ii(un')dfl' = / Ij(un)27rsin0d0 


ftu 


TT 

0 D (u) = J Ii(un')dfT = J~ Ij(ufl)27rsin0d0 


k/2 


Then the average upward and downward radiances 
are <I>u = 4>u/2tt, and <1 > d = 0 d / 277. Integrating 
equation (2.1 1) over fi U( the upward hemisphere of fi, 
and replacing l by its appropriate average gives 

1 d$u(u) w 

- — = ~4>\}{U) + — (0u(U) + 0 D (u)] + 

2 du 2 

w 

j-e u/ ^ + 2y 2 B(u) (2.12a) 
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Doing the same over fl D gives 


or 


1 d</> D (u) w 

-- -^ 7 — =-0 D (U)+-[0u(U) + <^d(U)1 + 

2 du 2 

■©- 

0 

+ 

NJ 1 — 

o. | n. 

c 1 -©- 

0 

II 

0 

(2.16c) 

w 

1 — e u/ ^ + 2 v 2 B(u) 

(2.12b) 

The solution of equation (2. 1 5) with B(u) of the form 

2 


of equation (2.8) and these boundary conditions is 

The solution of equations (2.12) is facilitated by 



letting 


<t>(u) = A|e 27U + A 2 e u/W + A 3 + A 4 e u/L 

(2.17) 

<f>( u) = </>u(u) + <t> d(u) 

(2.13a) 

where 


A0(u) = <t> u(u) - </> D (u) 

(2.13b) 



or 


1 2/i 0 

A, = [Jw(l -f 2 n 0 ) - - - + 4 B c 

! + 



1 +7 4/i 0 2 7 2 - 1 


0u — [0 + A</)]/2 

(2.13c) 

2L 7 2 

(2.18a) 


4B,(1+2L) 1 

4L 7 - 1 

</> d =W>-A<*>]/2 

(2.13d) 

4 Mo 2 

(2.18b) 

Adding (2.12a) and (2.12b) gives 


A2=JW 4 2 2 , 
4mo 7 ~ 1 

1 dAc p 

- = - (1 - w)$ + jwe u/,Wo + 4y 2 S 

2 du 

(2.14a) 

> 

t>> 

II 

.&• 

CD 

O 

(2.18c) 


4L 2 7 2 

(2. 1 8d) 

and subtracting gives 


A 4 — 4B) , , 

4L 2 7 2 -1 


1 d <*> 

= - A <f> 

2 du 


Differentiating (2.14b) gives 


(2.14b) 


Now, in the two-stream approximation, / 4jr lidO' = </>. 
Using equation (2.10) and inserting equation (2.15) 
into equation (2.7) gives 


1 d2< ft_ dA( ^ 

2 du 2 du 




l(i,e,g) = — j |Jw[p(g) 


1 d<£ 

1 ]e u/Mo + <t> ~~ Ti 
4 du 2 


e u/ M°du 


Substituting this result into (2.14a) gives 

1 d 2 <f> _ 7 

~ - 7 2< ^ + Jwe u/Mo + 4 7 2 B (2. 1 5) 

The boundary conditions on <p are that the radiance 
vanishes at u = -°o 


<£{-°°) — 0 (2.16a) 

and that there are no diffuse sources above the surface, 
so that 


Inserting equations (2.17) and (2.18) into this last 
equation and integrating gives 


l(i,e,g) : 


47TJU 


Jw(p -!)(— + -)-' + 
Mo M 


1 


A,W(2y + -)' 1 + A 2 


4mo 2 - 1 1 1 . 


W 


(— + -)-' + 
Mo M 


4L 2 - I 1 I ■' 
A 3M + A 4 — — (7 + - ) 


4L 


L n 


4>d ( 0 ) = I0(O)-A0(O)]/2 = 0 


(2.16b) 


which, after a modest amount of algebra, can be put 
into the form 
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(i.e.g) = J — M °~ [p(g) + H(w,/i 0 )H(w,/Li) -l] + 

4tt Mo + M 


Thus, the final form lor the radiance received by the 
detector is 


B 0 B| L 

(— 7H(w,m) + 7 2 H(w,L)H(w,m)1 (2. 1 9) 

7 r 7T L +/i 


where 


H(w,x) 


1 + 2 x 
1 + 2 yx 


( 2 . 20 ) 


•(i.e.g) = J 


W Mo 


47T Mo f M 


P(g)0 +' 


B c B | L 

yH(w,M) h — 

n ir L+m 


+ (l/k)tan(g/ 2 ) 


-) + 


H(W,Mo)H(W,m) - 1 


7 2 H(w,L)H(w,m) 


(2.24) 


More rigorously, the H-function defined by equa- 
tion ( 2 . 20 ) is actually an analytic approximation to a 
function defined by the integral equation 


w 

H(w,x) = 1 + — xH(w,x) 
2 


/ 


H(w,y) 
x +y 


dy 


( 2 . 21 ) 


The values of this function are tabulated in Chandra- 
sekhar (1 960). The differences between the H function 
defined by equation ( 2 . 21 ) and the approximation 
(2.20) are in all cases less than 4%, which is adequate 
for most remote-sensing applications. Also, it can 
readily be shown that 


/ 


H(w,x)dx 


1 + y 


( 2 . 22 ) 


One further phenomenon must be considered when 
the medium is illuminated by a highly collimated 
source: the opposition effect. This effect is a sharp peak 
in the brightness of the scattered radiance at zero 
phase angle due to the hiding of extinction shadows 
by the particles that cast the shadows. It occurs only 
when the particles are large compared to the wave- 
length, and acts only on the singly scattered part of 
the radiance. Its derivation is beyond the scope of this 
chapter, but is discussed in detail in Hapke ( 1 986). The 
result is that p(g) must be multiplied by a term of the 
form 1 + b(g), where 


Equation (2.24) describes the directional radiance 
of a particulate medium illuminated by a collimated 
source and with a thermal emission function that can 
be approximated by a constant term plus an expo- 
nentially decreasing term. This radiance is the quantity 
that would be measured by a field spectrometer or 
by a detector on a spacecraft. In equation (2.24), the 
first term on the righthand side describes the scattered 
radiance received by the detector, the term propor- 
tional to p(g) describes the singly scattered light, and 
the remainder describe the multiply scattered light. 
The second term on the right describes the thermally 
emitted radiance. 

However, the quantity that is frequently measured 
in the laboratory is the hemispherically integrated flux, 
which is the radiant power emitted per unit area into 
the entire upward direction. This quantity is given by 


n/2 


i 


$(0 - I i(u = o.aiMdn 


■/ 


l(0,e)M27rsinede 


A useful approximation for <f>(i) for the case of isotropic 
scatterers, p(g) = 1 , may be derived as follows. Be- 
cause the opposition effect is so narrow, the power 
included in it is negligible compared with the total 
power radiated into the entire hemisphere and may be 
ignored in calculating 4>(i). Thus, using the two-stream 
approximation and equations (2.13c) and (2.16b), 


it/ 2 


4>(i) — d>u( 0 ) J cose sine de = </> u( 0)/2 = c/>( 0)/2 


b(g) = 


1 + (l/k)tan(g/ 2 ) 


(2.23) 


b 0 is the amplitude of the opposition effect and k is 
its angular half-width. These two parameters depend 
on the nature of p(g), porosity, and particle size dis- 
tribution, and are discussed extensively in Hapke 
(1986). 


Substituting equations (2.17) and (2.18) into the last 
expression and performing a bit of algebra gives 


4>(i) = Jmo 


1 - y 


■ + 


27 


B 0 + 


2y 2yL 


1 + 27m 0 1+7 I + 7 1 + 27 L 


B, (2.25) 


4>(i) is the quantity that would be measured by an 
integrating sphere in the laboratory. It is also one of 
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the factors that determines the temperature of a 
surface in radiative equilibrium. The first term on the 
righthand side of equation (2.25) is the power per unit 
area scattered into the entire upward hemisphere by 
the medium; the second and third terms are the ther- 
mal power emitted into this hemisphere. 


2.2. REFLECTANCES, EMISSIVITIES, 
AND KIRCHHOFFS LAW 

2.2.1. Bidirectional Reflectance and BRDF 


2.2.3. Hemispherical-Directional Reflectance 

It will be shown below that the directional reflec- 
tance of a medium of isotropic scatterers illuminated 
by a diffuse radiance I 0 uniformly distributed over a 
hemisphere is often of interest in the measurement 
of thermal emissivity. Since the incident illumination 
is not collimated, the opposition effect can be ignored. 
The emergent radiance I e (e) under these conditions 
may be calculated from the reflectance portion of 
equation (2.24), with p(g) = 1 and b 0 = 0, by letting 
J = I 0 dfl = l 0 27rsinidi and integrating over i 


Several other quantities of interest for the interpre- 
tation of laboratory and remotely sensed data can be 
derived from equations (2.24) and (2.25). They will be 
discussed in this section. 

The coefficient of the incident illuminance J in 
equation (2.24) is the bidirectional reflectance 


. W Mo 

r d d(i-e,g) = 

4?r m 0 + M 


{[1 + b(g)]p(g) + H(w./i 0 )H(w.m) - 1 } 


(2.26) 


where b(g) is given by equation (2.23) and H(w,x) by 
equation (2.20). The reflectance is controlled primar- 
ily by w, the average single scattering albedo of the 
medium. 

The bidirectional reflectance distribution function 
(BRDF) is the ratio of the radiance scattered from a 
surface to the radiant power Jm 0 incident on a unit 
area of the surface. Hence, the BRDF is 


w 1 

BRDF(i,e,g) =— - — — 
47 t m 0 + M 


{[1 + b(g)]p(g) + H(w, Mo )H(w,M) - 1 } (2.27) 


n/2 


i 


w Mo 


I e (e) = I I 0 H(w,m 0 )H(w,m) 2 7r sini di = 


4ir Mo + M 


w r M 

— I 0 H(w,m) I H(w,m 0 )[1 ]dMo ; 

2 J Mo + M 


1 

W C 

I 0 {~ H(w,m) I H(W,Mo)dMo 
0 

I 


w r H(w,^ 0 ) 

- [- mH(w,m) j dMol } (2.29) 

2 J Mo +M 

o 


Using equations (2.21) and (2.22) gives 


I e (e) = I 0 {(1 - 7 )H(w,m) - [H(w,m) - 1] } 
loll -yH(w,M)] 


The hemispherical-directional reflectance is the re- 
flected radiance per unit incident radiance 


r hd (e) = I e (e)/I 0 = 1 - 7«(w,m) = , ~~~ (2.30) 

1 + 2yn 


2.2.2. Directional-Hemispherical Reflectance 

The directional-hemispherical reflectance is the 
fraction of power incident from a specific direction 
scattered by unit surface area into all directions into 
the upward hemisphere. The power per unit surface 
area of collimated illuminance J incident from a di- 
rection making an angle i with the normal is Jm 0 Hence, 
the directional-hemispherical reflectance is the co- 
efficient of Jmo in equation (2.25) 


r dh(i) 


1 ~7 
1 + 27Mo 


(2.28) 


where equation (2.20) has been used for the H func- 
tion. Note that r nd has the same functional depend- 
ence on e as r dh does on i. 

2.2.4. Bihemispherical Reflectance 

The bihemispherical reflectance is a theoretical 
quantity that is probably unobservable in practice but 
will be shown below to be important in the meas- 
urement and interpretation of thermal emittance. An 
approximate expression for the bihemispherical re- 
flectance can be found by a two-stream solution to 
the equation of radiative transfer in the form of equa- 
tion (2.12) for isotropic scatterers. Since we are inter- 
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ested only in reflectance in this section, B = 0. It is 
assumed that no collimated sources are present, so 
J = 0. Diffuse radiant flux </> 0 is assumed to be incident 
on the medium from the entire hemisphere above the 
surface, so that boundary condition (2. 1 6b) is replaced 
by c/> D (0) — <t>o- By the very straightforward calculation, 
the two-stream approximation for the bihemispherical 
reflectance is found to be 

1 -y 

r hh = 4>u(0)/4> o = (2.31) 

1 + y 

The bihemispherical reflectance is mathematically 
simple, yet is sufficiently accurate that it can be used 
for quantitative estimates in many radiative transfer 
problems. In support of the last statement, note from 
equation (2.28) that r dh (60°) = r hh , and from equa- 
tion (2.30), r hd (60°) = r hh . Also the bidirectional reflec- 
tance of a Lambert surface is, by definition, r L (i) = n 0 /tr; 
hence, if the particles scatter isotropically (p(g) = 1 ) 
and g is large enough that the opposition effect is 
negligible, then r dd (60 o ,60°,g)/r L (60°) = r hh . 

2.2.5. Directional Emissivity 

When the temperature T of the medium is constant, 
B, =0, and the radiance from a blackbody at that 
temperature is given by B 0 (T)/7r. Hence, from the 
thermal emittance part of equation (2.24), the direc- 
tional emissivity of the medium is 

1 +2u 

e d (e) = yH(w, M ) = y — — — (2.32) 

l + 2y^ 

If a large temperature gradient is present the emit- 
tance of the medium is given by the last two terms 
of equation (2.22). The blackbody emittance of the sur- 
face is then B[T(0)1 = (B 0 + If L= EA > 1, then 

py^yH(w,L)~ 1, and the emittance is given approx- 
imately by € d (e)(B 0 +B,)/rr. In this case the emittance 
corresponds to that of a blackbody at the temperature 
of the surface multiplied by the directional emissivity. 
As L decreases, the effects of B, decrease and the 
emittance is influenced more and more by the emis- 
sion corresponding to a body at the subsurface tem- 
perature B 0 . However, since A can hardly be smaller 
than the thickness of a monolayer of particles, which 
is of the order of 1/E, L can never be significantly 
smaller than unity. 

2.2.6. Hemispherical Emissivity 

When the temperature of the medium is constant, 
so that B| = 0, the hemispherically integrated flux 


from a blackbody at that temperature is B 0 . Hence, 
from the emittance part of equation (2.25) the hemi- 
spherical emissivity is 

2y 

e h = - (2.33) 

1 + y 

If there is a subsurface temperature gradient, then 
as with the directional emittance, if L> 1, the hemi- 
spherical emittance is « h (B 0 -I- B,) = e h B[T(0)|, and the 
effects of B 0 increase as L decreases. 

2.2.7. Kirchhoff's Law 

Kirchhoff's Law states that the sum of the reflec- 
tance and emissivity of a surface must be unity. It is 
of interest to inquire as to the sense in which Kirch- 
hoff's Law holds in a particulate medium. First, from 
equation (2.4f), 

6 = 1 -w (2.34) 

This expression, which shows that the sum of the aver- 
age particle emissivity and single scattering albedo 
equals unity, is simply a statement of Kirchhoff's Law 
for individual particles in the medium. 

Because the reflectance is often more convenient 
to measure in the laboratory than the emissivity, the 
reflectance is frequently measured and the emissivity 
calculated assuming Kirchhoff's Law is valid. How- 
ever, because a particulate surface can be character- 
ized by several different types of reflectances and 
emissivities, it is not always obvious which are the 
appropriate quantities to use. 

Comparing equations (2.30) and (2.32) it is seen that 

e d (e) = 1 - r hd (e) (2.35) 

Hence, in principle, if it desired to measure the direc- 
tional emissivity at some angle e, the surface should 
be flooded with diffuse light from all directions and 
the hemispherical-directional reflectance measured 
at the same angle e. However, in practice, it is usu- 
ally easier to measure the directional-hemispherical 
reflectance. Because r hd (e) and r dh (i) have the same 
dependence on their respective arguments, the direc- 
tional-hemispherical reflectance r dh (i = e) can be used 
in equation (2.35) to obtain e d (e). 

Similarly, comparing equations (2.31) and (2.33) 

t h = 2y/(l +y)= I - (I - y)/(l + y) = 1 - r hh (2.36) 

Thus, to measure the hemispherical emissivity the 
bihemispherical reflectance must be measured. This 
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is difficult, if not impossible, to do in practice. How- 
ever, as discussed in section 2.2.4, an equivalent value 
can be obtained by measuring r hd , r dh , or r dd /r L at the 
appropriate angles. 

2.3. DISCUSSION 
2.3.1. Emissivity 

In this section the properties of several of the quan- 
tities derived in this chapter will be discussed, along 
with some applications of interest in remote sensing. 
In order to avoid complications it will be assumed here 
that the particles scatter isotropically. The parameter 
that controls both the scattering and emissivity prop- 
erties of the surface is the single scattering albedo w. 
This parameter is a function of particle size, compo- 
sition, and wavelength. 

The variation of t h with t, or equivalently, of r hh with 
w, is shown in Fig. 2.3. Note that t h is a monotonic, 
nonlinear function of t. The slope of the curve is very 
large for small e, but flattens for e near 1 . Thus, a high 
spectral contrast is observed when e < 1 (w ~ 1 ), but 
the contrast is much smaller when e ~ 1 (w 1 ), a 
point emphasized by Conel (1969). Unfortunately, be- 
cause of the presence of restrahlen bands, most plan- 
etary regoliths have low albedos and high emissivities 
in the thermal infrared, so that the spectral contrast 
that can be observed is limited. 



e 


Fig. 2.3. Plot of hemispherical emissivity e h vs. particle 
emissivity « = y 2 and bihemispherical reflectance r hh vs. 
particle albedo w. 



Fig. 2.4. Plot of directional emissivity e d and hemispherical- 
directional reflectance r dh vs. angle of emergence e. 


The dependence of e d (e) and r hd (e) on emission angle 
e is plotted in Fig. 2.4. When e = 1 (w 1 ), e d — 1 at 
all angles, and the surface emits like a blackbody. 
Hence, for dark surfaces the assumption that the 
emissivity is independent of angle is a good approx- 
imation. As e decreases, the dependence of e d (e) on 
angle increases. 

2.3.2. Band Contrast in Reflectance 
and Emittance 

It is of interest to compare the amount of contrast 
in an absorption band when observed in reflectance 
and in emittance. From equation (2.26) it is seen that 
the bidirectional reflectance is linearly proportional to 
w(A), so that an absorption band with modest contrast 
in w(A.) will cause an equal or larger contrast in re- 
flectance. From equation (2.32) the directional emis- 
sivity is 


V 1 - w(\) (1 +2 m) 

e d (A,e) = — ===== 

1 + 2\J 1 - w(A) p 

The amount of contrast in t d depends on w and the 
particle size. If w(A) < 1, as is often the case, then the 
change in V 1 - w(A) with X will not be very large. In 
addition, this quantity appears in both the numerator 
and denominator of e d , which further reduces the con- 
trast. 

Reflectance and emittance from media in which the 
particles are smaller than the wavelength is poorly 
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understood. However, Salisbury and Easles ( 1 985) have 
demonstrated that the number of particles per unit 
volume can have a dramatic effect on both the mag- 
nitude of the reflectance (and thus on the emissivity) 
and on the spectral contrast within a band. A very 
porous powder of fine particles has low reflectance 
and spectral contrast. This can be understood if a small 
particle in a low-density powder scatters light approx- 
imately as if it were isolated. In that case, Q s cc(D/ 
A) 4 , while Q e °c (D/A) 1 , where D is the diameter {Kerker, 
1 969). Hence, w = Q s /Q e ^ (DA) 3 1 . Thus, e d (A,e) — I 

and the medium radiates practically like a blackbody. 

As another example, the effect of thermal emission 
on the contrast in absorption bands in the mid-infrared 
will be calculated. Suppose the spectrum of an asteroid 
surface at 3 AU whose regolith is at a temperature 
of 210 K is measured at e = 0°. The asteroid is 
illuminated by sunlight with an effective temperature 
of 5770 K at angle i = 45°. Assume that w = 0.40 at 
the center of an absorption band and w = 0.50 in the 
adjacent continuum. 

Using the equations developed in this paper it is 
readily calculated that, if observed only in bidirectional 
reflectance, the band-center-to-continuum ratio is 
0.72. If the band is at A = 4 pirn, the thermal emittance 
contribution to the radiance increases this ratio by 0.04 
to 0.76, so that the contrast is 2446, and the band is 
easily observable. However, if the band is at A = 5 pm, 
the thermal emittance increases this ratio by 0.24 to 
0.96, so that the band contrast is only 4%, and the 
band is difficult to observe. 


2.3.3. Radiative Equilibrium Temperature 


It is almost always assumed that the emitted radiance 
is independent of e, the angle at which the surface 

is observed. 

A more exact calculation using the equations of 
this paper is as follows. The fraction of visible sun- 
light reflected is the directional-hemispherical albe- 
do, so that the amount of sunlight absorbed per unit 
area of surface is (from equation (2.28) Jp 0 (l -r dh ) = 

Jm o TvH(w v ,p 0 ). 

To simplify the derivation, it will be assumed that 
the temperature of the medium is constant. However, 
it should be kept in mind that this is an approximation; 
a more rigorous calculation would take account of the 
-fact that the actual temperature increases into the 
medium Then from equation (2.33) the amount of 
energy radiated into the upper hemisphere is 2B| R y, R / 
(I + y IR ), where B, R is the Stefan-Boltzmann function 
aT 4 and T is the actual surface temperature. 

Equating the energy absorbed to that radiated gives 

B,r = oT 4 = )vMoy v H(w„,Mo)( 1 + 7iR)/2yi R 
so that the radiative equilibrium temperature is 


T = 


JvMo7vH(W v ,Mq)(1 + Tir ) 1 
2ay m 


1/4 


(2.38) 


If this surface is observed at angle e the 1R radiance 
is 


Iir(I.S) — c d (e)B| R — 7 V ( 

rr 2 tt 


+ 7v)Mo h (W v ,Mo) h ( w ir.M) 


(2.39) 


Finally, an expression for the radiative equilibrium 
temperature of a particulate medium with a smooth 
surface will be derived. In what follows, the subscript 
v will denote a quantity integrated over the spectrum 
of visible sunlight, and the subscript 1R will denote a 
quantity integrated over the thermal infrared. 

The usual expression for radiative equilibrium is 
found by equating the amount of visible sunlight ab- 
sorbed, taken as J v p 0 (l - A v ), where J v is the illumi- 
nance of visible sunlight and A v is a visual albedo, to 
the amount of infrared energy radiated, taken as e 
i R ctT c 4 , where e, R is an infrared emissivity, a is the 
Stefan-Boltzmann constant, and T c is the calculated 
surface temperature. This gives 


T c = 


f JvMo ( 1 A v ) i 
CTe IR 


1/4 


(2.37) 
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